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摘要 : 嫦娥 四 号 着 陆 器 将 搭载 低频 射电 频谱 仪 在 月 球 背面 进行 低频 射电 天 文 观测 。 该 
氏 频 射电 频谱 仪 的 观测 波段 为 0.1MHz~40MHz。 根 据 着 陆 器 在 中 国 空间 技术 研究 院 的 微波 
室 进 行 的 电磁 兼容 性 (Electromagnetic Compatibility ,简称 EMC ) 试验 结果 ， 着 陆 器 平 
台 在 该 频段 内 自身 存在 非常 强 的 噪声 ， 其 强度 甚至 小 没 掉 大 部 分 来 自 太 阳 爆 发 的 信号 ， 难 以 
探测 有 效 信 号 ， 实 现 预 期 的 科学 目标 探测 。 本 文通 过 模拟 仿真 分 析 谱 减法 ， 自 适应 滤波 以 及 
维 纳 滤波 这 三 种 方法 对 着 陆 器 噪声 消除 的 效果 ,从 而 选择 更 为 有 效 的 噪声 消除 方法 ,为 低频 
射电 频谱 仪 在 轨 探 测 任务 的 数据 处 理 提供 依据 。 
关键 词 : 低频 射电 频谱 仪 ; 嫦娥 四 号 ; 谱 减 法 ; 维 纳 滤波 ; 自 适应 滤波 
中 图 分 类 号 : P161 


0 引言 


美国 的 著名 工程 师 卡 尔 。 央 斯 基 在 1930 年 用 他 自己 建造 的 “旋转 木马 ”首次 观测 到 了 


来 自 银河 系 的 射电 辐射 "2， 自 此 展开 了 对 射电 天 文 的 广泛 研究 。 经 过 天 文学 家 的 研究 发 现 ， 


不 同 的 天 体 因 为 辐射 机 制 的 不 同 会 产生 不 同 波段 的 电磁 辐射 。 因 此 , 对 宇宙 进行 全 波段 观测 


一 件 非常 有 意义 的 寻 


旦 是 ， 地 球 大 气 的 存在 使 得 很 多 射电 波段 的 信号 特别 是 频率 低 于 


pum 
o 


10MHz 的 其 低频 波段 难 


a 


穿 透 到 地 面 。 然 而 ， 低 频 波段 的 探测 对 于 研究 天 体 的 演化 与 起 源 有 


着 至 关 重 要 的 作用 。 为 了 弥补 甚 低频 波段 的 观测 空白 ,国际 上 许多 国家 或 组 织 发 射 过 卫星 到 
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空间 进行 观测 ， 比 如 WIND 和 Ulysses. 


由 于 地 球 以 及 众多 在 轨 卫 星 的 电磁 干扰 , 空间 低频 射电 观测 也 面临 挑战 。 因 为 月 球 自 身 


LI 


可 以 遮挡 来 自 地 球 方向 的 射电 干扰 ， 所 以 月 球 背 面 的 射电 宁静 区 是 十 分 理想 的 观测 地 点 。 欧 


航 局 CESA) F 2015 年 1 月 提出 了 在 2025 年 对 月 球 背 面 发 射 一 颗 载 有 低频 射电 频谱 仪 的 着 


陆 器 的 FARSIDE 计划 中 。 我 国 提出 嫦娥 四 号 探 月 计划 。 嫦 娥 四 号 是 世界 首 颗 在 月 球 背 面 软 


着 陆 和 近视 探测 的 航天 器 。2018 年 5 H 21 H, WRU S REE “RSE” BOA, Abus 


和 巡视 器 将 于 年 底 发 射 ,该 着 陆 器 上 搭载 了 一 台 科 学 载荷 低频 射电 频谱 仪 , 是 人 类 在 月 球 背 


面 首次 实现 对 宇宙 和 太阳 的 低频 段 〈0.1MHz~40MHz) 进行 观测 。 


根据 着 陆 器 整 器 在 微波 上 暗室 的 EMC 试验 结果 ， 昌 然 来 自 外 界 的 信号 干扰 可 以 被 很 好 的 


屏蔽 掉 , 但 是 着 陆 器 本 身 搭载 的 电子 设备 以 及 其 它 的 载荷 对 低频 观测 存在 电磁 干扰 , 由 于 无 


法 对 平台 进行 工程 上 的 调整 , 所 以 只 能 通过 数据 处 理 的 方式 来 消除 平台 的 干扰 。 目 前 在 语音 


信号 以 及 雷达 信号 处 理 等 方面 ， 降 低 信 号 的 噪声 处 理 主要 有 三 种 方法 : 谱 减 法 、 自 适应 滤波 


和 维 纳 滤波 方法 。 本 文 将 结合 低频 射电 频谱 仪 的 科学 目标 ， 对 比分 析 信 号 增强 效果 ， 得 到 最 


优化 的 方法 ， 为 低频 射电 频谱 仪 在 轨 探 测 任务 的 数据 处 理 提供 依据 。 


1 低频 射电 频谱 仪 的 工作 原理 以 及 噪声 来 源 


低频 射电 频谱 仪 系统 由 四 根 接收 天 线 〈 三 根 Sm 长 的 长 天 线 A、B、C, 一 根 20cm 长 的 


短 天 线 D)， 前 置 放 大 器 以 及 电子 学 单元 组 成 。 三 根 长 天 线 A、B、C 是 相互 垂直 的 5 Chu 


下 图 1 所 示 )， 接 收 空间 的 电磁 信号 的 三 个 方向 上 的 分 量 ， 根 据 电磁 的 传播 理论 就 可 以 计算 


出 接收 的 信号 的 强度 以 及 方向 。 短 天 线 D 与 长 天 线 B. C 处 于 同一 个 平面 ， 根 据 电磁 波 的 


传播 理论 以 及 不 同 长 度 的 天 线 的 增益 的 不 同 ，A、B、C 长 天 线 可 以 同时 接收 到 来 自 宇宙 空 
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间 的 远 场 信号 和 来 自 着 陆 器 本 身 的 近 场 信号 ， 而 短 天 线 D 因 天 线 短 且 离 着 陆 器 近 ， 因 此 收 


到 主要 是 来 自 近 场 的 着 陆 器 噪声 信号 口 。 


图 1 低频 射电 频谱 仪 四 根 天 线 吕 


Fig.1 The four antennas of the very low frequency radio spectrometer 


ERU ^ E P 300 BS IUE E A SCF VE RS £9 8 AIK 3E (LOOK Hz~2MHz) Al re 251 


[3] 


(1MHz~40MHz) 两 个 频段 ， 频 谱 仪 的 主要 技术 参数 如 下 表 1 所 示 ， 
表 1 低频 射电 频谱 仪 主要 技术 参数 门 


Tab.1 Main technical parameters of Very low frequency radio spectrometer"! 


3i H 设计 参数 实测 结果 
工作 频率 100KHz~40MHz 100KHz~40MHz 
接收 机 灵敏 度 nV/VHz) 6 8.1 
低频 : 5KHz 4.88 KHz 
频率 分 辨 率 
高 频 ，100KHz 97.66 KHz 
最 大 数据 率 (bps) 4.5M 4.34M 
功 耗 (W) 24 24 


在 低频 射电 频谱 仪 工作 的 同时 , 着 陆 器 上 其 他 部 分 载荷 以 及 一 些 电子 模块 也 同时 处 于 运 


行 状态 ， 比 如 电源 控制 器 、 数 传 调制 器 、 发 射 段 RE102、 数 传 控 制 SMU 等 。 下 图 2 给 出 了 
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月 球 表面 的 电磁 波 通 量 密度 , 而 天 线 附 近 区 域 的 电场 强度 和 接收 机 输出 的 电压 值 可 以 根据 下 


面 公式 计算 得 出 ; 


G G i IL 
_ “preamp "Receiver coax 
Uout = 14+Zq/Z, *he E ( 1 ) 
a/4preamp 


AP, Uon 为 接收 机 的 输出 电压 值 ，Gpreamp 和 Greceiver 分 别 是 前 置 放 大 器 增益 和 接收 机 增益 ， 


ILeon 为 高 频 电费 插 损 ,Ze 和 Zpreamp 分 别 是 天 线 阻抗 和 前 置 放大 器 阻抗 ,为 天 线 有 效 长 度 ， 


巨 则 是 天 线 处 的 电场 强度 。 
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图 2 月 球 表面 的 电磁 波 通 量 密度 站 


Fig.2Radio emission levels on the surface of moon 


公式 〈1) 中 的 各 项 参数 可 以 根据 低频 射电 频谱 仪 生产 方 提 供 的 定 标 报告 获 取 ， 结 合 


[2] 


项 参数 以 及 公式 〈1) 可 以 将 接收 机 测 得 电压 数据 反 演 出 空间 电场 以 及 通 量 密度 。 图 2 中 给 


出 了 各 天 体 到 达 月 球 表面 的 通 量 密度 , 图 中 的 着 陆 器 噪声 是 根据 在 中 国 空间 技术 研究 院 的 微 


波 暗 室 当 中 对 嫦娥 四 号 着 陆 器 整 器 进行 的 EMC 试验 中 测 得 工 况 三 (着 陆 器 电源 控制 嚣 主 份 、 


数 管 管理 单元 主 份 、 数 传 调制 器 主 份 、 载 荷 电 控 箱 主 份 、 低 频 射电 频谱 仪 等 依次 加 电 连 续 工 


VE) 天 线 A 接收 到 的 高 频段 数据 反 演 得 到 。 从 图 中 可 以 发 现 着 陆 器 的 噪声 通 量 密度 大 于 普 


通 的 太阳 爆发 产生 的 通 量 密度 , 与 太阳 爆发 峰值 通 量 密度 相当 。 所 以 不 对 接收 到 信和 号 进行 降 


噪 处 理 ， 我 们 难以 反 演 出 太阳 爆发 的 真实 强度 。 
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2 低频 射电 频谱 仪 有 效 信号 提取 方法 


由 于 低频 射电 频谱 仪 的 三 根 主要 天 线 A、B、C 收 到 的 目标 信号 里 都 会 包含 有 来 自 着 陆 


器 自身 的 噪声 信号 的 干扰 , 所 以 在 对 信和 号 进行 提取 使 用 之 前 需要 通过 降 品 算法 来 抑制 噪声 并 


最 大 限度 的 分 离 出 有 用 的 信号 。 用 x(n) 表 示 观 测 到 信号 ，s(n) 表 示 有 用 信号 ，v(n) 表 示 品 声 


信和 号。 它们 中 间 的 关系 可 以 表示 为 : 
x(n)=s(n)+v(n) (2) 


在 本 文 的 模拟 仿真 中 ， 将 一 组 强度 与 太阳 爆发 强度 接近 的 模拟 信号 作为 有 用 信号 s(n)， 


用 EMC 试验 测 得 的 工 况 三 的 着 陆 器 噪声 当 作 噪声 信号 v(n)， 通 过 自 适 应 滤波 ， 谱 减法 以 及 


维 纳 滤波 的 方法 来 进行 仿真 处 理 。 


2. 1 谱 减 法 


谱 减法 是 信号 增强 的 一 个 常见 方法 , 其 原理 是 在 频 域 上 用 信和 号 加 噪声 的 幅 值 减 去 噪声 的 


幅 值得 到 有 用 信和 号 的 幅 值 。S.EBoll 于 1979 年 改进 了 该 方法 ， 称 为 改进 的 谱 减 法 口 。 改 进 的 


谱 减 法 在 谱 减 的 过 程 中 使 用 的 是 噪声 的 平均 功率 ， 可 以 有 效 得 减 小 剩余 的 “音乐 噪声 ”。 改 


进 谱 减 法 可 以 用 公式 表示 为 


|IX(n)|” —a*D(n); Xm >a* Din) 
icr = | b*D(n) |X(n)l" «a*D(n) (3) 


其 中 Da) 为 噪声 的 平均 幅 值 ， 


D(n) = 元 Zi1DiCDI(4) 


公式 (3)、(4) 中 Y 等 于 1 或 2， 为 1 时 就 是 幅 值 相 减 ，y 为 2 时 就 是 功率 相 减 。a (a>1) 


为 过 减 因子 , bb>0) 是 补偿 因子 ， 可 以 通过 调节 两 个 参数 来 减少 周期 性 的 误差 。 结 合 低频 射 


电 频 谱 仪 的 结构 特点 以 及 对 科学 数据 的 要 求 ， 可 以 用 长 天 线 A、B、C 收 到 信号 的 幅 值 的 平 
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均值 减 去 短 天 线 D 收 到 信号 的 幅 值 的 平均 值 来 获取 信号 的 幅 值 。 


由 于 在 设计 上 天 线 A、B、C 与 天 线 D 的 长 度 以 及 材质 不 一 样 导 致 其 增益 也 不 相同 ， 所 


以 在 进行 谱 减 的 时 候 需 要 补充 不 同 增益 带 来 的 信号 强度 差异 。 具体 来 说 , 先 取 一 小 段 时 间 内 


的 信号 计算 出 它们 的 对 消 系数 。 然 后 在 后 续 的 处 理 过 程 中 将 D. 天 线 收 到 信和 号 的 平均 幅 值 乘 


以 对 消 系数 来 获取 噪声 幅 值 。 可 以 结合 低频 射电 频谱 仪 的 结构 特点 来 运用 改进 的 谱 减 法 进行 


降 噪 处理 ， 以 天 线 A 和 天 线 D 为 例 可 以 用 公式 表示 为 


MAORORAA EA 
Scot - do IPC) * Vo Ol; VACA < IPC * vo PY 
其 中 VA(D 是 长 天 线 A 收 到 的 信号 进行 储 里 叶 变 换 后 在 相应 频 点 的 平均 值 , Vp(D 则 是 短 天 线 


D 收 到 的 信号 在 相应 频 点 的 平均 值 , y 取 1。P(f) 通 过 取 太 阳 宁 静 时 的 一 小 段 时 间 内 四 个 接收 


机 接收 到 的 信号 分 别 算出 的 ， 


_ Val) 
P,(f) = V (6) 


pb 人 


处 理 流程 如 图 3 所 示 。 


D 天 线 收 到 信号 做 
FFT 并 计算 平均 幅 
fi 


RF fä v (n) 


取 一 段 时 间 内 D 天 
线 收 到 信号 做 FFT 
并 计算 平均 幅 值 


A、B、C 天 线 的 
幅 值 减 去 D 天 线 
WAE RELAXT 
系数 后 的 值得 
到 信号 幅 值 


EINHEDR Hi 
对 消 系数 P 


取 同 一 段 时 间 内 


平均 幅 值 


Hepa fos E 


一 一 带 品 信号 x (n) 


图 3 改进 谱 减 法 运算 示意 图 


Fig.3Improved spectral subtraction diagram 


下 面 三 种 数据 处 理 方法 以 天 线 A 与 天 线 D 数据 进行 仿真 为 例 ， 谱 减法 仿真 结果 如 图 4 


所 示 ; 


100 
一 一 带 噪 信号 x(n) 
谱 减 法 处 理 后 信号 
80 一 一 有 用 信号 s(n) 
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图 4 改进 谱 减 法 仿真 结果 


Fig.4Simulation results by improved spectral subtraction 


e 图 4 中 展示 的 是 FET 后 的 频 域 信号 ， 横 轴 表 示 频 率 〈 单 位 为 MHz) , 纵 轴 为 幅度 〈 单 位 


= 为 dBuV)。 图 4 中 黑色 的 曲线 是 带 有 噪声 的 信号 ， 蓝 色 的 曲线 是 经 过 改进 谱 减 法 处 理 后 的 


信号 ,红色 的 曲线 是 有 用 信和 号。 从 图 4 中 可 以 看 出 该 方法 在 整个 波段 对 噪声 都 有 有 效 的 抑制 ， 


在 10MHz 以 下 的 频段 内 噪声 的 幅度 下 降 了 20dB 左右 , 在 其 它 频段 也 能 将 噪声 降低 10dB Æ 


A 右 ， 且 处 理 后 的 信号 波形 基本 能 还 原 原 本 的 有 用 信号 ， 在 背景 噪声 较 低 的 频段 


(20MHz~40MHz) 处 理 效 果 更 好 。 


2. 2 维 纳 滤波 


AAR D 中 可 知 ， 我 们 进行 信号 处 理 的 目标 是 得 到 不 含 噪声 的 有 用 信号 s(n). {AE 


在 实际 的 信号 处 理 过 程 中 ， 我 们 求 得 的 有 用 信号 并 不 完全 等 于 s(n)， 只 能 是 sq gr ex 


估计 值 , 因此 对 信号 的 处 理 可 以 看 作 是 对 s(n) 的 估计 , 处 理 信 号 就 是 找到 一 个 最 佳 的 估计 器 。 


维 纳 滤波 就 是 一 个 寻找 最 佳 估计 器 来 估计 信号 的 方法 ”1。 


y(n) 
s (n) +v (n) 


图 5 维 纳 滤波 器 
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Fig.5 Wiener filter 


图 5 为 维 纳 滤波 的 原理 示意 图 ， 其 中 h(n) 是 维 纳 滤波 器 的 滤波 参数 ，y(n) 是 对 有 用 信号 


s(n) 的 估计 值 。 用 公式 可 以 表示 为 


y(n) = x(n) * h(n)(7) 


根据 最 小 误差 准则 AMS) fie = E[{s(n) — y? p^ $3 ^it n] EAE HB de EEE BL, 


Els(n) * x(n)] 


hn) = 一 一 一 一 一 一 (9) 
E[x’(n)] 
将 上 式 取 傅 里 时 变换 可 以 导出 ， 
P. (k 
H k — SX 
(k) P W (9) 


其 中 PKA s(n) 与 x(n) 的 互 功 率 谱 密度 ，Px(k) 为 x(n) 的 功率 谱 密 度 。 由 于 有 用 信号 s) 
噪声 信号 v(m) 互 不 相关 ， 即 有 Pw)=0, 则 可 以 得 到 
Pox (k) = P.(k) (10) 


P; (k) = P;(k) + Pa(k) (11) 


将 式 (10) 和 式 C11) WAR (9) 中 可 以 得 到 式 (12) 


1 
Hk) = 1 - — (12 
(k) aad 


AUKO2) B, y CO) AERE AE oh DML rti i I RR x t a f E iT 


根据 天 线 D 和 天 线 A 接收 到 数据 可 以 计算 出 y(k)， 从 而 计算 出 最 佳 的 滤波 参数 。 下 面 图 6 


为 仿真 结果 。 


作 期 刊 


电压 /dbuV 


一 一 带 噪 信号 x(m) 
一 一 维 纳 滤波 处 理 后 的 信号 
一 一 有 用 信号 s(m) 


qr Wye 


从 图 6 中 可 以 看 出 维 纳 滤波 也 可 以 降低 噪声 的 影 


1 1 1 Å 
20 25 30 35 40 
频率 /MHz 


L 
15 


图 6 维 纳 滤波 仿真 结果 


Fig.6Simulation results byWiener filter 


响 ,但 是 相对 于 谱 减 法 的 降 品 效 果 要 差 


Ab E 


些 。 


(0-20MHz) 仍 难以 分 辩 


2.3 自 适应 滤波 降 品 


基本 也 能 还 原 原 有 用 信号 的 波形 ， 但 是 处 理 后 的 信号 在 


噪声 比较 大 的 波段 


Fi 


出 有 用 信 


< 不 同 于 维 纳 滤波 , 自 适应 滤波 器 可 以 根据 收 到 信号 不 断 地 自动 的 调节 自身 的 参数 来 最 大 
= 限度 的 降低 噪声 信号 中。 自 适应 滤波 已 经 广泛 运用 在 系统 识别 、 信 道 均衡 、 信 号 增强 以 及 信 
号 预测 方面 。 本 文中 主要 运用 的 是 它 的 信号 增强 的 功能 , 它 的 结构 示意 图 如 下 图 7 所 示 。 自 
适应 滤波 算法 至 少 需 要 两 套 接收 设备 ,其 中 一 套用 来 接收 带 有 噪声 的 信号 , 一 套 专门 用 来 接 

收 噪声 ， 然 后 通过 自 适应 算法 求 得 有 用 信号 ， 而 低频 射电 频谱 仪 正好 有 这 两 套 接收 设备 。 


x(n) 


vi(n) e(n) " 


É | 适 UN ye 波 At "m 


适应 滤波 器 结构 图 


图 7 


Fig.7Structure of the adaptive e filter 


如 上 图 7 所 示 为 一 个 自 


适应 降 品 滤波 器 的 结构 示意 图 ， 


中 包括 有 ) 


的 信号 s(n) 以 及 噪 
号 , e(n) 表 示 误 差 信 和 号。 它们 之 间 的 关系 可 以 表示 为 


e(n) = 


式 (13) 中 w(n) 表 示 的 是 滤波 器 的 加 权 系 数 ， 该 系数 可 以 通过 


其 


声 信 号 v(n))， 而 vin) Wet A 


中 x(n) 是 指 接收 机 收 到 信号 ( 


N 


一 个 接收 机 接收 到 的 噪声 信 


x(n) — XL [wv (n — {)]3) 


最 小 均 方 误差 准则 (LMS) 或 者 最 


小 二 乘法 RLS) 评估 Ele*(n)] 来 调整 ， 使 滤波 器 工作 在 最 佳 的 状态 n ”J。 下 面 利用 最 小 二 
乘法 (RLS) 来 进行 仿真 。 
100 T 
ARS x(n) 


一 一 自 适应 滤波 处 理 后 的 信号 
一 一 有 用 信号 s(n) 


] | M 
Ea mw 


Fig.8Simulation results byadaptive filter 


从 图 8 中 可 以 发 现 利用 自 适 应 滤波 的 方法 仿真 的 结果 比较 差 ， 降 品 效 果 不 明 显 ,， 对 信和 号 


的 还 原 度 也 不 高 。 在 噪声 信号 中 一 些 幅 


适应 滤波 处 理 


) AA 


后 噪声 信号 幅度 没有 下 降 。 


3 低频 射电 频谱 仪 背 景 噪声 消除 结果 分 析 
2.1 到 2.3 已 经 通过 三 种 不 同 的 数据 处 理 方法 对 着 陆 器 噪 


度 比较 大 的 频 点 处 (如 1.73MHz、2.32MHz、2.9MHz 


声 进 行 了 仿真 处 理 


Es 
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应 滤波 方法 是 在 时 域 的 处 理 方法 , 改进 谱 减 法 和 维 纳 滤 波 则 是 在 频 域 的 处 理 方 法 。 表 1 给 出 


了 这 三 种 信号 增强 办 法 对 信 


言 号 信 噪 比 的 提升 效果 的 对 比 。 


de 2 三 种 方法 处 理 数据 前 后 的 信 品 比 


Tab.2 SNR before and after the data processing by three different algorithms 


数据 处 理 方法 


处 理 前 信 噪 比 /dB 


处 理 后 信 噪 比 /dB 


=> 


售 噪 比 提升 /dB 


改进 谱 减 法 维 纳 滤波 9 适应 滤波 (RLS) 
-27.64 -27.64 -27.64 
1.99 -4.64 -25.35 
29.63 23.00 2.29 


从 表 Y 中 我 们 可 以 发 现 改 进 谱 减 法 和 维 纳 滤波 的 处 理 结果 明显 好 于 自 适 应 滤波 , 且 其 中 


改进 谱 减 法 的 信 喉 比 提 升 是 略 大 于 维 纳 滤波 的 。 目 适应 滤波 信号 增强 系统 降 噪 能力 受 参 考 输 


入 端的 噪声 信号 vi 与 带 品 信 号 中 的 噪声 信号 的 相关 度 有 很 大 的 影响 。 在 低频 射电 频谱 仪 


设计 的 时 候 用 于 接收 信号 的 长 天 线 A、B、C 与 用 于 接收 着 陆 器 噪声 信号 的 短 天 线 D 的 各 项 


参数 不 同 , 以 及 仿真 的 信号 中 有 用 信和 号 的 信 噪 比 非常 的 低 等 因素 导致 了 自 适 应 滤波 降 噪 的 效 


果 变 差 。 长 天 线 A. B. C 与 短 天 线 D 的 参数 不 同 导致 D 天 线 与 长 天 线 收 到 噪声 信号 的 不 同 


可 以 在 频 域 上 通过 评估 短 天 线 与 长 天 线 在 不 同 频 点 处 的 增益 差异 来 抵消 , 所 以 在 频 域 处 理 的 


两 种 算法 在 降 噪 时 能 取得 比较 好 的 效果 。 


为 了 验证 低频 射电 频谱 仪 的 性 能 以 及 数据 处 理 方法 , 我 们 在 中 国 科 学 院 国家 天 文 台 密云 


观测 站 食堂 楼 顶 安装 好 低频 射电 频谱 仪 以 及 接收 天 线 进行 了 外 场 试验 , 并 且 在 距离 接收 天 线 


100 米 左 右 的 宿舍 楼 里 用 牺 


F 意 波形 发 生 器 连接 天 线 发 射 模拟 信号 。 而 着 陆 器 噪声 信号 是 采用 


的 EMC 试验 工 况 三 的 试验 结果 ， 由 一 个 放 在 接收 天 线 下 方 的 1 米 长 的 小 天 线 连 接 波形 发 生 


器 发 出 。 下 图 9 为 其 中 一 组 试验 天 线 A 接收 到 的 数据 在 利用 改进 的 谱 减法 处 理 后 的 结果 ， 


其 中 发 射 的 有 用 信号 为 1.8MHz 的 点 频 信和 号 。 
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图 9 改进 谱 减 法 处 理 结果 


Fig.9Noise reductionby improved spectral subtraction 


图 9 中 截取 了 1.7MHz-1.9MHz 的 信号 ， 从 图 中 可 以 看 出 经 过 改进 谱 减 法 处 理 后 ， 着 陆 


器 的 噪声 被 有 效 的 抑制 了 ， 幅 度 下 降 了 10dB， 而 有 用 信号 的 幅度 并 没有 被 降低 ， 信 品 比 也 


从 -25.23dB 提升 到 了 -2.60dB 。 通 过 改进 谱 减 法 处 理 后 的 信号 ， 信 噪 比 提升 了 -22.63dB ， 说 明 


改进 谱 减 法 在 外 场 实验 中 对 于 着 陆 器 背景 噪声 信号 的 降低 也 是 有 效 的 。 


在 月 球 背 面 进行 极 低 频 (0.1~40MHz) 探测 是 非常 有 科学 意义 的 1, 同时 由 于 该 频段 信 


号 微弱 信 品 比 低 ， 所 以 提取 有 用 信号 也 是 十 分 困难 的 。 本 文 利用 EMC 测试 的 着 陆 器 噪声 以 


及 嫦娥 四 号 搭载 的 低频 射电 频谱 仪 对 月 球 背面 能 接收 到 的 目标 信号 强度 进行 了 仿真 实验 。 利 


用 改进 的 谱 减 法 、 维 纳 滤波 以 及 自 适应 滤波 这 三 种 常见 的 信号 增强 的 算法 对 仿真 信号 进行 了 


提取 ,其 中 改进 的 谱 减 法 对 信号 的 提取 效果 显著 , 经 处 理 后 信 品 比 有 明显 的 提升 ,可 以 直观 


的 从 处 理 后 的 波形 中 看 见 目 标 信号 波形 。 并 且 改 进 的 谱 减法 对 在 密云 进行 的 外 场 试 验 的 数据 


TA 


的 处 理 也 有 很 好 的 效果 。 因 此 , 可 以 将 改进 的 谱 减 法 作为 低频 射电 频谱 仪 主要 的 预 处 理 方法 ， 


实现 对 太阳 爆发 进行 较为 准确 的 观测 。 低 频 射电 频谱 仪 在 月 球 背面 进行 探测 的 过 程 中 , 温度 


等 环境 的 变化 可 以 引起 对 消 系数 的 改变 , 所 以 在 实际 探测 过 程 中 需要 结合 实际 环境 情况 , 每 
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Research of The Method forReducing Background of Very 
Low Frequency Radio Spectrumon Chang’E-4 
Zhang Tao" ^ ?, Su Yan" * 3 
(1. University of Chinese Academy of Sciences, Beijing 100049, China, Email: zhangtao@nao.cas.cn; 
2. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China; 
3. Key Laboratory for Exploration in the Moon and deep Space, Chinese Academy of Sciences, Beijing 100101,China) 
Abstract: Chang’E-4 will carry a Very Low Frequency Radio Spectrometer(VLFRS), which 
will land on the far-side of the moon forRadio observation. The VLFRS has an observation 


frequency of 0.1] MHz-40MHz.According to Electromagnetic Compatibility(EMC) test results in 


the Microwave darkroom of China Academy of Space Technology, the lander platform is quite 


201811.00086v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


noisy in this band, and the intensity of the noise is even stronger than that of most signals from the 


sun burst. In this pape, we will simulate and optimize background elimination through improved 


spectral subtraction, adaptive filtering and wiener filtering, and choose a more effective method 


foreliminating background. 


Key word: Very low frequency radio spectrometer;Chang’E-4;Improved spectral 


subtraction; Wiener filtering;Adaptive filtering 


